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ABSTRACT 

We present near-infrared VLT ISAAC imaging and spectroscopy of the peculiar Virgo galaxy 
NGC 4438, whose nucleus has been classified as a LINER. The data are supplemented by 
mid-infrared imaging, and compared to previous WFPC2 HST broadband images. Images 
and position-velocity maps of the [Fe ll] and H2 line emissions are presented and compared 
with the distribution of the optical narrow-line region and radio features. Our results show that 
shocks (possibly driven by a radio jet) contribute to an important fraction of the excitation of 
[Fell], while X-ray heating from a central AGN may be responsible for the H2 excitation. 
We address the question whether the outflow has an AGN or a starburst origin by providing 
new estimates of the central star formation rate and the kinetic energy associated with the 
gas. By fitting a Sersic bulge, an exponential disc and a compact nuclear source to the light 
distribution, we decomposed NGC4438's light distribution and found an unresolved nuclear 
source at 0.8" resolution with Mk — —18.7 and J — H = 0.69. Our measured bulge velocity 
dispersion, 142 km s^^, together with the standard A^bh — f relation, suggests a central black 
hole mass of log(Albh/M0) ^ 7.0. The stellar kinematics measured from the near-infrared 
CO lines shows a strong peak in the velocity dispersion of cto ^^178 km s^^ in the central 
0.5", which is possible kinematic evidence of a central black hole. We calculated a general 
expression for the integrated Sersic profile flux density in elliptical geometry, including the 
case of 'disky' isophotes. 
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1 INTRODUCTION 

NGC 4438 is a large peculiar spiral galaxy, with a disturbed stel- 
lar disc and an even more heavily disturbed interstellar medium 
(ISM). Located near the centre of the Virgo cluster, NGC 4438 has 
undergone a violent collision with the nearby giant elliptical M86 
( [Kermey et al. 12008 1, and may also be experiencing ongoing ram 
pressure stripping due to an interaction with the Virgo intraclus- 
ter medium ( [Vollmer et al.|[2009j l. NGC4438's nucleus has been 
classified as a particularly interesting LINER (low-ionisation nu- 
clear emission-line region, |Heckman| 1980) 1. The spectra of LINERs 
are characterized by the presence of emission lines from atomic 
species of low ionisation state. By definition, LINERs are galax- 
ies which host nuclei with emission-line ratios that satisfy the fol- 
lowing criteria: [O III] A5007/H/3 < 3, [O I] A6300/Ha > 0.05, 
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and [N II] A6583/Hq > 0.5 ( |Osterbrock & Fer land'200yi. Their 
emission line spectra are similar to those observed in narrow-line 
regions (NLRs) of gas in Seyfert 2 galaxies. However, some LIN- 
ERs have relatively powerful central black holes, and tend to have 
broader emission lines (e.g., NGC 1052, |Ho et al.|1997| l. 

The most likely mechanism to explain the excitation in these 
objects is photoionization either from an active galactic nucleus 
(AGN) or from a strong stellar continuum jHo et al.|2003| . AGN 
photoionisation models fit both the low-ionisation spectra of LIN- 
ERs and the high-ionisation spectra of the Seyfert NLRs with sim- 
ilar nuclear emission, but different nebular conditions (such as the 
electron density and the incident ionising luminosity; see discus- 
sion by Ho et al._2003,,psterbrock & Ferland 2005) . 

The nucleus of NGC 4438 lies at the root of a nuclear bub- 
ble, expanding to the north-west (NW), which has been imaged 
in Ha+[N II] emission by Kenney & Yale|p002 ). In contrast with 
other similar systems, it is not clear what is powering this bub- 
ble, since there is neither a strong starburst nor a strong AGN in 
NGC 4438. For example, the galaxy NGC 2782 posseses a central 
starburst (extending over a radius ~200 pc) which provides the me- 
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Figure 1. A possible decomposition of Ha A6562 and [N II] AA6548, 6583 
emission, using tlie minimum number of components (a linear baseline, 3 
narrow lines and 1 broad component - or 13 free parameters). The x-axis is 
wavelength in A, while the ?/-axis is relative intensity. The solid circles are 
the data points, the solid line is the best fit model, and the dotted lines are 
the contributions from each component. Data courtesy of L. Ho. 



chanical luminosity that drives its central winds fjogee et al.'1999'). 
The Seyfert galaxy NGC 3079 harbours both a powerful AGN and 
a circumnuclear starburst ( [Veilleux et al.|1994l >, while M82's out- 
flows are driven only by a nuclear starburst ( [Lehner et al.|1999| l. 
Moreover, there are many other objects powered by composite sys- 
tems where a circumnuclear starburst coexists with a central low- 
luminosity AGN ( [Rodriguez-Ardila et al.|[2b05[ l. [Levenson et al.| 
([2003) observed the Seyfert 2 starburst galaxy NGC 5135, distin- 
guishing both the AGN (unresolved) and the starburst (spatially ex- 
tended over ~200 pc). 

Broad Ha emission in NGC 4438 has been tentatively inferred 
by |Ho et al^ ( |1997| l from the fitting of optical spectra. We have 
repeated the fitting of the Ha emission line profile on their data 
but with a smaller number of Gaussian components than |Ho et| 
[aT] ( [1997j >, obtaining similar results. Fig. [T| shows that a broad, 
FWHM~2050 krns~^, Gaussian component is clearly present 
in the emission-line complex. This spectral feature is thought to 
be indicative of the presence of active nuclei harboring a broad 
emission-line region ( |Ho et aL|i997) ). 

Near-IR (near-infrared) light can escape high opacity and 
dusty environments more easily than Ha photons can, making it 
possible to detect heavily obscured line-emitting regions. Near-IR 
emission-lines; such as [Fe II] 1.257 ^m and H2 2.122 /im, are less 
sensitive to extinction than optical lines such as Ha. Near-IR obser- 
vations can provide better extinction estimates, and can yield new 
information on the intensity of a potentially obscured central AGN 
or starburst. 

Results of a near-lR spectroscopic survey of LINER galax- 
ies carried out by |Larkin et al.|p998^ , showed that the [Fell] line 
{a*D-ji2 — 0^-^9/2 at 1.257 /im) is the most commonly detected 
emission-line in LINERs (see also, [Rodriguez^ Ardila et al.|2005l l. 
The H2 line (1 — 0S(1) at 2.122 /im) is also a common feature. 
On the other hand, the Pa/3 emission-line (or H I 3-5 at 1.282 /im), 
which is the strongest near-IR recombination line available at low 
redshift in AGN and starburst galaxies, is only found in emission in 
20% of the LINER galaxies in |Larkin et"ar] ( |1998| l sample. Br7 (or 
H I 4-7 at 2. 166 /im) is undetected in all the galaxies of the survey. 



Table 1. ISAAC imaging observing log and calibration factors. 



Date 


Filter 


Integration 


Zero point 




SeeinJ^ 






(sec) 


(mag) 


(mag) 


(arcsec) 


broad-band data 


Apr 14 2003 


J 


400 


24.85 


0.025 


0.86 


Apr 14 2003 


H 


400 


24.31 


0.016 


0.62 


Jan 13 2003 




284 


24.12 


0.010 


0.58 


narrow-band data 


Apr 14 2003 


1.26 


600 


21.88 




0.69 


Apr 14 2003 


1.28 


600 


21.69 




0.56 


Jan 13 2003 


2.07 


800 


21.62 




0.55 


Jan 14 2003 


2.13 


800 


21.80 




0.57 


Jan 13 2003 


2.17 


800 


21.76 




0.51 



Foreground extinction (see Section 3.1 . 

Effective seeing at the observed wavelength, given by the FWHM of the 
best-fitted Gaussian profile. 

NGC4438's combination of peculiar spectral features, com- 
plex morphology and the possibility of studying the nuclear bub- 
ble at different wavelengths led us to obtain new near-IR imag- 
ing and spectroscopic data. In this work, we present the results of 
ISAAC imaging and spectroscopy of NGC 4438, focusing on anal- 
ysis of line emission maps and on two-dimensional modelling of 
the galaxy surface brightness. We want to study the morphology of 
NGC 4438 and whether there is a nuclear point source embedded 
in the bulge. In Section [2] we begin by giving a brief description 
of the observations. In Section[3]we discuss the nuclear line emis- 
sion based on both imaging and spectroscopy, together with ex- 
tinction estimates. In Section [33] we derived the line-of-sight stel- 
lar velocity dispersion from CO bands absorption features. In Sec- 
tion[4]we present the 2D decomposition of the surface brightness of 
NGC 4438. In Section[5]we discuss and compare our findings with 
previous studies in the X-ray, optical and radio wavelengths, focus- 
ing on the energetics of the nuclear source and the surface bright- 
ness modelling results. Appendix [A] gives a detailed description of 
the data reduction of the imaging data. In AppendixjBjwe describe 
the surface brightness model used to fit the light distribution and 
how we measure structural parameters. 

Throughout this paper we use a distance to NGC 4438, 
near the centre of the Virgo cluster, of 16 Mpc (at which 
1" corresponds to a distance scale of 77.6 pc). All the 
data reduction and analysis were carried out using the Perl 
Data Language (http://pdl.perl.org, [Glazebrook et| 
[aT] [T997J. The stellar kinematics was derived using IDL 
(http : / / www . ittvis . com/idl/). 



2 OBSERVATIONS 
2.1 Imaging 

We observed NGC 4438 with the VLT ISAAC imager and spec- 
trograph (Moorwood et al.[|1998| l, which has a field of view of 
152" X 152" and a pixel scale 0.148" pixePV Our data con- 
sist of the following broad and narrow band (NB) filters: J, H, 
Ks, NB1.26 ^m, NB1.28 fim, NB2.07 ^im and NB2.13 ^m filters. 
Our observing strategy involved acquiring 3 to 10 s exposures in 
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Table 2. Spectroscopy observation log. 



Date 


Filter 


Spectral domain 

{p,m) 


Slit width 
(arcsec) 




Integration 
(sec) 


(degrees) 


Notes on quality and spectral features 


Apr 04 2002 


SL 


2.80-4.00 


0.6 


600 


24x0.94 


-57.10 


Good. No PAHs. 


Apr 04 2002 


SK 


2.25-2.37 


0.3 


8900 


2x300 


+107.25 


Good. CO band-heads 


Feb 09 2003 


SK 


2.11-2.23 


0.3 


8900 


2x300 


+107.25 


Good. CO band-heads 


Feb 09 2003 


SK 


2.11-2.23 


0.6 


4400 


2x300 


-57.65 


Bad, fringing 


Apr 14 2003 


J 


1.25-1.30 


0.6 


5200 


2x600 


-57.65 


Regular. Resolved [Fe II] emission 


May 23 2003 


J 


1.25-1.30 


0.3 


10500 


2x300 


-37.55, +17.25, +76.65 


Good. Resolved [Fe II] emission for PA-37.55 


May 23 2003 


SK 


2.11-2.23 


0.3 


8900 


2x500 


-37.55, +17.25, +76.65 


Regular. Resolved H2 emission for PA-37.55 



" Spectral resolution from ISAAC manual. 

^ Position angle on the sky, convention is positive from North to East. 




arcsec 

Figure 2. Overlay of the ISAAC 120-arcsec slit on the J-band ISAAC im- 
age of NGC 4438 for position angles: +17.65, -37.65, -57.65. Convention is 
positive position angle from North to East. The x- and j/-axes ai'e east and 
north offsets in arcsecs from the nucleus at 12000 RA 12'^ 27™ 45.6= and 
J2000 DEC +13° 00' 32" {from HST astrometry). 



2.1.1 Complementary imaging data 

We also observed NGC 4438 witli the TIMMI2 camera on the 
ESO 3.6m telescope on UTC 2002-01-01X08:50, with a lens scale 
of 0.3 arcsec per pixel. The target was acquired in the NB 11. 9 ^m 
filter by blind pointing on the HST coordinates of the nucleus, with 
a 5" rms positional uncertainty in a field of view of 90" x 60". The 
object is detected in NB 11.9 /im, which is free of emission lines. 
Yet our aligned chopping and nodding strategy, with a throw of 30" 
North-South, fits one of the negative images on the array. The pres- 
ence of the chopped image confirms our detection. We enhanced 
the signal-to-noise ratio of our image by smoothing with a Gaus- 
sian kernel, with a dispersion of 2 pixels, or 0.6". Flux calibration 
was obtained by comparison with HD81797, with a 11.9 /im flux 
density of 100.57 Jy. By fitting an elliptical Gaussian to the stan- 
dard star we infer an angular resolution of 0.76 x 0.71" FHWM. 
After smoothing, the resolution of the image in Fig. [lO] is about 
0.96". We set the astrometry by tying the centroid of the 11.9 pim 
emission to the HST coordinates for the nucleus. 

We complemented the infrared data with recalibrated WFPC2 
HST images, broadbands F450W (B), F675 (i?) and F814W (J). 
We retrieved these data from the HST archive located at the Cana- 
dian Astrophysics Data Centre (CADC). The CADC pipeline recal- 
ibrates the images with up-to-date calibration files. The details of 
these observations are described by |Kenney & Yalel p002l. Since 
the pixels covering the nuclear emission were saturated in R and / 
images, we used only the B-band image to carry out the photome- 
try of the nucleus. 



a dithered pattern consisting of four frames with vertical and hori- 
zontal offsets of ~ 22" and ~ 9", respectively. 

Data reduction and calibration procedures are summarized in 
Appendix |A] Table [T| gives the dates, integration times and cali- 
bration parameters for each observation. Fig. |2] shows the ISAAC 
J-band image of NGC 4438 after flat-fielding, bad-pixel correction, 
and sky and bias subtraction. 

Only the central region of the field (inner 40 arcseconds) has 
a symmetric point-spread function (PSF). The outskirts of the ar- 
ray show extended tails due to optical aberrations. In this central 
region we report a seeing of ~ 0.5" for the NB filters, while for the 
broadband filters we report a seeing of 0.6" — 0.8" (see Table[T](. 



2.2 Medium-resolution near-IR spectroscopy 

Near-IR spectra of NGC 4438 were obtained with ISAAC at VLT 
on 2003 March-April. We acquired medium resolution spectra at 
central wavelengths of 1.274 /im, 2.170 /^m, and 2.310 /im, and 
slit widths of 0.3 and 0.6 arcsec. The list of observations, filters, 
position angles, spectral domains, resolutions and slit widths can 
be found in Table|2] Fig.|2]shows some of the slit positions overlaid 
on the J-band ISAAC image. 

The sky background (including OH skylines) was removed by 
differencing along the slit, with nod throws of 50 arcsec (or 30 arc- 
sec in the case of standard stars, STDs hereafter). Wavelength cal- 
ibration was obtained by comparison with arc lamps. The spec- 
tra were extracted in 0.74 arcsec centred on the peak of emission. 
Telluric absorption spectra and flux density calibration were ob- 
tained by observing the early-type star HD 115709 (spectral type 
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Figure 3. ISAAC near-IR spectra of the nucleus of NGC4438. Top spec- 
tram is centred on 1.274 /xm, while the bottom spectrum is centred on 
2.170 /im. Top and bottom spectra were acquired with position angles are 
-37.55° and -57.65°, respectively. Notice the absence of the hydrogen re- 
combination lines, Pa/3 at 1.28 and Br7 at 2.16 ^m, in the spectra. 
Both spectra were extracted from a 0.74 arcsec apperture. Wavelengths are 
in microns, flux densities are in units of 6.5xlO~^^ W fim~^ . 




arcsec 

Figure 4. Extinction in the central 40" of NGC4438, superposed on 
HST Ha+[NII] contour map jKenney &^Me 20021. The highest ex- 
tinction coiTesponds to the region surrounding the nucleus, with a peak 
of Av ~ 4.4 mag at ai'ound l"east of the nucleus. Contours are - 
0.01, 0.06, 0.27, 0.62, 1.1, 1.73, 2.49, 3.4, 4.44, 5.62, 6.94 times 2.4 x 
10~^® ergs~^ cm~^. The x- and y-axes are east and north offsets in 
arcsecs from the nucleus at J2000 RA 12^ 27™ 45.6= and J2000 DEC 
-f 13° 00' 32" (from the HST astrometry). 

the collapsed 0.6 arcsec slit is 3 X lO"^'^ W m"2/.im"^ at 2.8 /im, 
and 6 x 10^^° W ^m~^ at 4.0 /xm. The noise varies as a func- 
tion of wavelength, between 10^^^ to 10"^^ W ^mT^ . This 
spectrum shows no evidence of the presence of polycyclic aromatic 
hydrocarbons (PAHs) emission in NGC4438. 



AlIV) and comparing with a black-body spectrum. Stellar absorp- 
tion features at Pa/3 and Br7 were accounted for by fitting Voigt 
profiles in a total wavelength range of 0.015 /im in the vicinity of 
the stellar absorption lines. Samples of the reduced spectra, centred 
at 1.274 and 2.170 /im, are shown in Fig. [3] 

Flat-fielding was hampered by short-term variations in the 
CCD pixel gains. The observatory pipeline flat fields could not sup- 
press the low-level features on scales of 0.005 seen in the spec- 
tra of Fig. [3] The features under ~ 2 x 10^^^ W m~^ /im~^ were 
not reproduced in different exposures, which suggests flat-fielding 
artifacts. Try as we might we could not improve on the pipeline 
flats. We treat the low-level artifacts as noise, so that the depth of 
the spectroscopy did not meet our expectations. The ISAAC near- 
IR spectroscopy is nonetheless informative on the kinematics of the 
CO band-heads, [Fell] 1.26 /im and H2 2.12 /im. 



2.3 Low-resolution spectroscopy 

We acquired a low-resoution L-band (3.8 /^m) spectrum at PA 
-1-56.94, in 7 ABBA cycles, for a total integration time of 1 h (with 
individual DITs of 0.935 s and 4 NDITs, a chopthrow of 20 arcsec, 
and a chopping frequency of 0. 1 1 Hz). The spectrum shows a con- 
tinuum rising towards longer wavelengths. This L-band spectrum 
does not show any spatially resolved features. The flux density in 



3 RESULTS 

3.1 Extinction estimates 

A consequence of the dependence of interstellar extinction on 
wavelength is that the near-IR imaging often reveals surprising 
morphological differences when compared to optical observations 
(see [Cardelli et al.|[T989| for a description of the Galactic ex- 
tinction law). Nuclear extinction estimates are very uncertain in 
NGC4438. The Balmer decrement Ha/H/3 gives = 3.4 (us- 
ing the extinction law from [Cardelli et al.|1989^ , but this represents 
a luminosity-weighted average over the central 2" x 4" extraction 
aperture of the low-resolution ground-based spectroscopy reported 
by |Ho et ar]jl997 [l, and also a lower limit to the true extinction. 

In order to further investigate the spatial distribution of the 
dust in the system we generated an extinction map (see Fig.[4|. This 
was done by transforming the J — K colour image into a colour ex- 
cess map, E{J — K). We assumed a flat and constant stellar pop- 
ulation, as expected for bulges of spirals, with E{J — K) = 0.8 
(calculated over an 7" x 7" aperture in a region devoid of structure). 
The colour excess map E{J — K) and the column of hydrogen nu- 
clei Nn are related by Nh/E{J-K) = 1.1 x 10^^ cm"^ mag"\ 
assuming a Galactic gas-to-dust ratio and Rv ~ 3.1 ( [Tokunag^ 
2000 1. We inferred a mean colour excess E( J — K) = 0.40 for an 
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aperture matched with the 2" x 4" used by Ho et al. ( 1997^ . Using 
the standard IR interstellar reddening law i Cardelli et al.|1989] ), it 
implies a visible colour excess of E{B — V)= 0.7, which is fully 
consistent with the colour excess reported by |Ho et al.| ( |199'7] l, of 
E{B — V) = 0.61. As revealed by the presence of strong dust 
lanes, the nuclear extinction is highly variable toward the central 
hundred parsecs ( [Kenney & Yale|2002[ l. 

The spatial distribution of extinction and dust is shown in 
Fig.|4] The lower edge of the Ha bubble emission is coincident with 
bands of strong extinction. These region of high extinction seem to 
be surrounding the nucleus (given by the peak of Ha emission), 
with a peak of Ay ~ 4.4 mag at around l"east of the nucleus. 

The measured colour excess implies an average hydrogen col- 
umn density TVh = 2.7 x 10^^ cm~^ at the base of the NW bubble 
(shown in Fig. |4|, and a mean column TVh = 1.2 x 10^^ cm~^ 
for the rest of the NW bubble (using the total to selective absorp- 
tion from ,Tokunaga 2000). This column density is in reasonable 
concordance with the values reported by Mac hacek et al.| ( (2004| ). 
|Machacek et al.| P004| l obtained a best-fitted column density of 
~ 2 x 10 cm~ from modelling of the X-ray spectra of the NW 
bubble. The uncertainties in our extinction estimates are mainly due 
to variations in the stellar population near the nucleus. 



3.2 Maps of extended emission 

Fig. [5] shows the [Fell] 1.26 /^m and H2 2.122 fj,m maps of line 
emission, after subtracting the continuum emission. The fact that 
there is some negative emission at the base of the NW bubble im- 
plies that there is some Pa/3 emission contaminating the continuum 
image. We are certain that this contamination does not affect the 
NW extended emission since the spectrum shows no trace of Pa/3 
emission at a noise level of 1 a = 1.06 x lO^^'' W mT^ /.im^^, 
when the slit is aligned with the NW bubble (PA -37.55 and PA 
-57.65, see Fig. [3j. In the case of H2 emission we used a scaled 
version of the broadband image to subtract the continuum. The 
presence of this strong [Fe II] emission is consistent with LINER 
spectroscopic surveys (Larkin et al.'1998 ), confirming the LINER 
classification of the NGC 4438's nuclear region. 

The outflow shells are very asymmetric in the optical lines and 
the radio and X-ray continuum. This provides further evidence that 
the outflow shells are intrinsically different, and that the observed 
dissimilarities are not due principally to extinction. In the near-IR 
we have detected only the NW outflow, with no trace of the south- 
eastern outflow neither in [Fe II] nor H2 emission. 



3.2.1 [Fe II 7 extended emission 

The emission-line map shown in Fig. [5^ indicates that [Fell] 
emission is coincident with HST Ha-l-[N II] emission-line de- 
tected by |Kenney & Yale| ([20021, Chandra X-ray emission re- 
ported by 'Machacek et al.'J2004^, and the complex radio contin- 
uum blob A described by Hummel & Saikia| ( |1991| ). Emission from 
a counter-shell in the southeastern (SE) region is clearly detected 
in Ha-l-[N II] and radio (see [ Kenney & Yale|[2002) , but is much 
fainter, and is located much further from the nucleus (9"vs 4") than 
the much-brighter NW shell. Fig. [5] also reveals a lack of this SE 
emission. This is because the SE emission is 15 times fainter than 
the NW component in HQf-l-[N II] ( [Kermey & Yale|2002[ l, and even 
weaker in the X-ray (a factor of 32 lower counts than the NW shell, 
[Machacek et al.|2004^ . Also, most of the bubble-like emission seen 
in Ha-l-[N II] corresponds to [N II] emission since Ha is about two 



times weaker than ionized nitrogen emission according to Fig. [T] 
The bright H II regions, detected in Ha-l-[N II] emission 8" toward 
the south of the nucleus by [Kenney & Yale|p002| ), do not appear 
neither in Pa/3 or [Fe II] emission nor in X-ray emission ( [Machacek] 
[et al.|2004| ). 

In AGN galaxies, the [Fe II] emission is thought to arise from 
the region radiating narrow line emission ( [Mouri et al.|2000][Rif-[ 
[fel et al. ||2006). Such regions in AGN can be produced either as 
a result of photoionization by a nuclear source, including X-ray 
heating from the central AGN, or via shock excitation by radio jets 
( [Rodriguez- Ardila et al.[[2005j [Riffel et al.[[2006 ). The dominant 
excitation mechanism of the [Fe II] emission is still under debate. 
[Simpson et al.[ jI996) have argued that photoionization is the dom- 
inant excitation mechanism of [Fe II] and that shocks by radio jets 
account for only about 20 per cent of the emission in Seyfert galax- 
ies. 'RodriguezjArfjlaet^(|2005^ have shown that the [Fe Il]/Pa/3 
ratio is a good indicator of the relative contribution of photoioniza- 
tion and shocks, since [Fe II] seems to be more tightly correlated 
with the radio emission than hydrogen recombination lines in AGN 
( [Simpson etal.|1996^ . 

Radio maps of NGC 4438 presented by [Hummel & Saikia[ 
( [1991^ show extended emission, with a shell-like morphology, that 
is well aligned with Ha and [Nil] (see figures 3 and 4 in [Kenney] 
[& Yale|2002[ l, and therefore also spatially coincident with [Fell]. 
This morphological correlation, along with the non detection of ex- 
tended Pa/3 emission in our spectra of the NW bubble, may imply 
that shock excitation by radio jets contributes to an important frac- 
tion of the [Fe II] emission in NGC 4438. 



3.2.2 H2 extended emission 

Molecular hydrogen emission is present in many host galaxies of 
AGNs ( [Rodrfguez-Ardila et al.[]2005] [Riffel et al.[|2008l ]2009|. 
In the traditional AGN picture, a dusty torus shields the H2 gas 
from the dissociating AGN radiation field. The H2 emission in 
NGC 4438 seen in Fig.jSjj could also be associated with the central 
AGN engine but it does not seem to trace the NLR gas ( [Rodriguez^ 
[Ardila et al.|2005| >. As it can be seen in Fig. ]5] the molecular hy- 
drogen is distributed only in the inner few hundred parsecs. The 
H2 emission is much stronger on the north side of the bubble than 
on the south side. Fig. [5|5 shows that most of the molecular hy- 
drogen gas is concentrated toward the circumnuclear region tracing 
the presence of a possible dusty torus, as in the case of NGC 3727 
( [Rodn'guez-Ardila et al.|]2005) or NGC 4051 ( [Riffel et al.[]2008l ). 
Also, some H2 emission extending towards the NW bubble can be 
seen in our observations. 

The position- velocity maps depicted in Fig. |6] show that the 
kinematics of the [Fe II] and the H2 emitting gas are distinct. Fur- 
thermore, the nuclear [Fe II] line is broader than the H2 emission. 
As measured from the nuclear spectra shown in Fig. [3] the full- 
width at half-maximum (FWHM) of the [Fell] and H2 A2.12^m 
lines are 415 and 270 km s~^, respectively. This implies that the 
[Fe II] is originating from a kinematically more disturbed gas than 
the H2 emitting gas, which is in agreement with previous ob- 
servations of other active galaxies, such as NGC 21 10 |Storchi-[ 
[Bergmann et al.|I999[ > and NGC 405 1 jRiffel et al.|2008) . A pos- 
sible interpretation for this is that the H2 emitting gas is more re- 
stricted to the galactic plane, perpendicular to the radio jets, while 
the [Fe II] emitting gas extends to higher latitudes from the galac- 
tic plane ( [Riffel et aL|2008[|Storchi-Bergmann et al.|1999l l. These 
results suggest that the H2 excitation is likely to be dominated by 
X-ray heating from the central AGN. 
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Figure 5. a) ISAAC [Fe II] continuum subtracted emission-line map of the central 10 " in NGC 4438. Note that there is a negative arc due to Pa/3 contamination 
in our continuum frame. The colour-scale is exponential and the contours correspond to Hq:-h[Nii] emission at 0.12, 0.14, 0.18, 0.27, 0.39 and 0.55 times 
32.9 MJy sr-^. b) Molecular hydrogen (H2) in the central 10" of NGC 4438. The colour-scale is exponential and the contours are 0.04, 0.09, 0.18, 0.31, 0.45 
and 0.62 times 21.9 MJy sr^^. The x- and j/-axes are east and north offsets in arcsecs from the nucleus. 
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Figure 6. Results of long-slit observations of NGC 4438. Panel a) shows [Fe II] at 1.257 /^m extended emission while panel b) corresponds to H2 at 2.122 /jm. 
Spectral PSF subtraction was carried out in order to remove the background and continuum of the galaxy. The slit position for both spectra coiTresponds to 
the direction along the major-axis of the NW bubble (position angle is -57.65 degrees), thus giving the velocity gradient across the shell. Intensity is given in 
units of the rms noise on a pseudo-colour exponential scale. The j;-axis is the spatial offset along the slit in arcsec, and the x-axis shows the velocity in the rest 
frame of the object in km s~ ^ . 



3.2.3 PaP extended emission 

Fig. [5^ also shows a negative arc to the east of the nucleus, perpen- 
dicular to the ISrW shell. In this arc Pa/3 is stronger than the [Fe II] 
emission. The distance between the nucleus and the arc is about 
0.5" (~35 pc), which is much larger than the values pointed out for 
the radius of the putative torus. These values are typically smaller 
than 5 pc (see, for example, |jaffe et al.|2004[|lVIinezaki et al.|2004| 
|Riffel et al.||2009] l. A possible interpretation for this emission is 
that it might be associated with circumnuclear star forination in the 
outer parts of the torus. 

The circumnuclear region seen in negative in Fig. [5^ has a Pa/3 
flux density of 0.5 mjy (integrated over the arc-shaped feature), 
which corresponds to a Pa/3 luminosity of 3.6 x 10*" erg s~ ^ . If this 



emission has a starburst origin and we assume solar abundances, 
and a Salpeter initial mass function i/i oc M~^ '^^ with an upper 
mass cutoff of 100 Mq, we can estimate the star formation rate 
(SFR) in the arc feature. |Kennicutt| |l998j find by extrapolation: 
SFR(]Vl0 yr"^) = 7.9 x 10"*^L(Ha) (ergs"^), computed for 
Case B recombination at Te = IC* K. Under these conditions the 
SFR as a function of Pa/3 luminosity is, 

SFR (]Vl0 yr"^) = 1.4 X 10""^ L (Pa/3) (ergs^^), (1) 

for jHa/ipa/3 = 17.6, whcrc j is the emissivity of the 
line ^Storey & Hummer||1995| l. The extinction-corrected nuclear 
star formation rate, using a central extinction Ay of 4 mag, is 
1 M© yr~^. It represents an overestimated SFR in the circum- 
nuclear region because this emission could include a contribution 
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Figure 7. Medium-resolution spectrum centred at 2.310 with a slits of 
0.74 arcsec. The red line represents the best galaxy model fitted by pPXF 
(see text) as a linear combination of template spectra. 
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Figure 8. Line-of-sight (LOS) stellar velocity dispersion (cr) as a function 
of distance along the minor-axis of the galaxy. The LOS velocity dispersion 
reaches a maximum at the nucleus of the galaxy with crmax ~ 180 kms-i 



the Gemini Near-IR spectral templates librarjj^ Stellar templates 
were rebinned to the spectral resolution of the galaxy spectra (15.47 
km pixel"^). The best-fitting parameters were determined by 
minimisation in the pixel space using the penalized pixel fit- 
ting method (pPXF) developed by |Cappellari & Emsellem| ( |2004| l, 
which required a signal-to-noise ratio of at least 30 in order to get 
reliable stellar kinematics. Errors were estimated by using a Monte- 
Carlo scheme. They were obtained as the standard deviation of 
the kinematical parameters (Vios and ctios) for many realisations 
(A'^ = 300) of the input spectra by adding Gaussian noise to the 
best model of the galaxy spectrum. The kinematical study of the 
nature of this putative central mass object requires detailed kine- 
matical model (e.g., [Cretton & van den Bosch|1999[|Cappellari et| 
|al.|2009|f . Unfortunately, the scarcity of our data does not allow us 
to perform this kind of analysis. 

The results shown in Fig.[8]that NGC 4438 exhibits a strong 
peak of the line-of-sight velocity dispersion along the minor axis 
on the centre of the galaxy. The velocity dispersion varies from 100 
to 180 km in a region of about 3"(~ 230 pc) in size, which 
could suggest the presence of a massive central object. 

The radius of the sphere of influence of this possible central 
massive object is just ~0.l", which is considerably smaller than 
our 0.74" spatial resolution. Therefore, the enhancement observed 
in the stellar velocity dispersion curve in the inner -^230 pc, where 
sigma rises from about 100 up to 177 km s~^, might be only due to 
the galaxy bulge potential. However, the sphere of influence is not 
a hard limit to the influence of a black hole on its host bulge. The 
influence of the compact object will still be felt at a larger radii but 
to a lesser extent. Comparing Fig. 8 with the dispersion profiles of 
other AGNs such as NGC 2549, whose sphere of influence is also 
~0.1" (see fig. 4 in |Krajnovic et al.|2 009), NGC4438's velocity 



dispersion seems to be even steeper in the central couple of arcsec. 
Therefore, we cannot rule out the presence of a supermassive black 
hole as responsible for the central enhancement in the velocity dis- 
persion profile of NGC 4438. 



4 SURFACE BRIGHTNESS MODELLING OF NGC 4438 



from the AGN, due to gas that has been ionised by an accreting 
black hole. For a global description of the star formation history of 
NGC 4438 as a whole see |Boselli et al.| ( (2005l >. 



3.3 CO-band stellar kinematics 

NGC 4438 is an inclined large spiral galaxy and its kinematics are 
undoubtedly difficult to study mainly due to the multiplicity of 
components that build up spiral galaxies. Moreover, NGC4438's 
highly disturbed morphology and the presence of dust make the in- 
ference of kinematic measurements a difficult task. On the positive 
side, its inclined disc permits a relatively simple identification of 
the direction which might be expected to define one of the princi- 
pal axes of the bulge. 

CO-band stellar line-of-sight kinematics along the minor axis 
(PA = -1-107.5°) were derived by fitting the galaxy spectra (see 
Fig.|7]( with a linear combination of template spectra chosen from 



4.1 Model description: general considerations 

Does the nucleus seen in our broadband images correspond to a 
cuspy surface-brightness profile, or is it a point source? Our ap- 
proach to addressing this question consists in modelling the stel- 
lar component by fitting the bulge with a r^^" Sersic law, and the 
outer region (disc) with an exponential profile. The nuclear com- 
ponent will be a priori represented by a point source at the cen- 
tre of the galaxy. The main advantage of this approach is that it 
lets us identify a nucleus over the stellar contribution, whereas a 
non-decomposition approach would try to reproduce all the profile 
without giving back information about separate components. The 
universality of the Sersic plus exponential disc decomposition al- 
lows us to compare with other objects. A detailed description of 
the surface brightness model, as well as the optimisation algorithm 
can be found in Appendix [B] 



^ see near-IR resources at: 

http : / /www . gemini . edu/sciops/ instruments/ 
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Table 3. Best NGC4438 galaxy structural parameters. 



Filter 


A 


q 


c 


a 


n 


Re 




Rh 




n' 


R'e 






log-Mbh 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 
















bull 


;e + disc + Sersic 


nucleus 












J 


30 


0.57 


1.10 


68.3 


1.65 


3.88 


71.76 


13.46 


42.65 


1.75 


~ 


572.63 


Tin 


7.19 


6.90±0.45 


H 


30 


0.57 


1.14 


68.4 


1.47 


3.53 


94.76 


12.77 


57.08 


1.93 


~ 


502.12 


21.8 


6.97 


6.68±0.45 


Ks 


30 


0.57 


1.13 


68.7 


1.52 


3.49 


73.39 


12.86 


41.34 


1.71 


~ 


570.92 


13.0 


6.66 


6.75±0.45 


B 


30 


0.50 


1.20 


68.9 


0.80 


2.50 


23.00 


8.70 


30.00 


1.70 


0.20 


154.00 


7.56 







Notes.- (1) Filter name. (2) Angular size of the fitted region. (3) & (4) generalised ellipticity parameters. (5) Inclination angle in degrees. (6), (7), (8), (9) and 
(10) are the Sersic structural parameters for the bulge, iJe and Rf^ are in units of arcseconds . (11 ), (12) and (13) are the Sersic structural parameters for the 
nuclear source. (14) Total bulge luminosity in units of 10^^ erg s calculated from Equation |B3| (15) Black hole masses obtained from the A^bh"luminosity 
density relation ^Marconi et al.|2003^. (16) Black hole masses obtained from the A4bh"™ relation ^Graham & Driver|2007j , the errors are estimated using the 
intrinsic uncertainties in the |Graham & Driver|j2007} relation. All radii are given in arcsecs, the intensity units are MJy sr"'^ and the black hole masses are 
expressed in Mq . 




4.2 Application to tiie near-IR ISAAC data 

We fitted the nuclear region (30" x 30") of NGC4438 with para- 
metric functions, as described in Appendix [B] The model is the 
sum of three components: a bulge, a disc and a compact source. 
The dust features present near the central region have been taken 
into account by correcting the images using the extinction map (see 
Fig.g. 

The best-fitting structural parameters for each filter are sum- 
marised in Table [3] The best fit is always obtained with n close 
to 1.7 and i?e ~ 3" for the bulge component present at near-IR 
wavelengths (see nomenclature in Appendix|Bj. The ellipticity and 
shape of the isophotes are quite constant with increasing radius, 
with an ellipticity e = 0.5 and slightly disky isophotes (c = 1.2, 
see Appendix [B|(. 

Since the integrated flux density presented in [Peng et al.| 
( |2002^ is undefined in the case of 'disky' isophotes, we calculated 
a general expression for the integrated Sersic profile flux density, 
in generalised elliptical geometry, including the case of 'disky' 
isophotes (Equation |B3]l. The integrated luminosities of the bulge 
were computed from the flux densities integrated over all radii. 

The result of the two-dimensional fit for the J— filter is pre- 
sented in Fig.|9] The left panel shows a cut along the major axis of 



the galaxy (dashed curve), the model (solid curve) and its compo- 
nents: bulge (dotted curve), disc (dash-dot-dot curve) and the nu- 
clear component (dash-dot curve). That figure shows an excellent 
agreement between the model and the galaxy's light distribution. In 
this central region we can see a conspicuous point-source standing 
out of the bulge. However, galaxies with power-law profiles may 
show substantial differences along the major and minor axis, due to 
isophote twists or ellipticity. For this reason, the right panel shows 
the surface brightness azimuthally-averaged along a radial elliptical 
annuluQ this panel also shows an excellent match between model 
and data at large radii. 

We find evidence for a nuclear point-source unresolved at 0.8" 
resolution (see Fig.|9| with Mk = -18.7 and J - AT = 0.69. Its 
extinction-corrected near-IR integrated fluxes in a 3" aperture for 
each band are listed in Table |4] The computed bulge luminosity in 
each broadband are listed in Table [3] A very similar system was 
studied by [Peng et aL] ( |2002^ , who found a nuclear point source 
embedded in NGC 4278, which is an elliptical Seyfert galaxy, also 
with a LINER nucleus. 



^ With the same nomenclature of the left panel, in logarithmic scale 
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Table 4. Near-IR nuclear source photometry. 



Filter name 


Central A 


Fx 




Lx 




(/xm) 


(mJy) 


(mag) 


(lO^'i ergs-i) 


J 


1.25 


10 


-17.96 


7.35 


H 


1.65 


12 


-18.66 


6.67 




2.16 


8 


-18.67 


3.40 



'^The specific luminosity was computed as Li, = iivD^uFi,, at a distance 
D = 16Mpc. 

4.3 Application to the i/Sr-WFPC2 data 

Since the HST B band image shows much more structure than the 
near-IR data (due to extinction being more severe at shorter wave- 
lengths), we followed two different approaches in order to perform 
the photometry of the nucleus]^ The first approach was similar to 
the one used for the near-IR data, i.e., fitting a Sersic bulge plus 
an exponential disc and a nuclear source. The nuclear source was 
represented by an extra Sersic profile, since the central source in 
the B-band image is clearly resolved. The result was a less promi- 
nent bulge in the optical than in the near-IR. In the optical case the 
surface brightness is well represented by a bulge with n = 0.8 and 
a compact source well-fitted by a steeper Sersic component with 
n = 1.7 and = 0.2", i.e,. a resolved nuclear source of approxi- 
mately 16 pc in size. An n value of 0.8 seems very low for a bulge 
of a spiral and it may not represent the real bulge surface distribu- 
tion, this is probably because of the severe obscuring effect of dust 
at those wavelengths. 

The second approach consisted of modelling the structure of 
the central 2" x 2" in the B-band image with a base of legen- 
dre polynomials Pi{x) and Pk{y), both at order l,k — 6, in or- 
der to get accurate photometry of the nucleus. We added a re- 
solved compact-source, represented by Fo 5{xcan, J/ccn) convolved 
with a gaussian of FWHM — 0.4". The free parameters were 
the (/ + 1) X {k + 1) legendre coefficients, the centroid of the 
compact-source (a^ccn, J/ccn), the FWHM of the nuclear source, and 
the central intensity Fo. The result was similar to the near-IR case: 
a resolved compact source with FWHM — 0.4", equivalent to 
31 pc, appears to stand out of the nuclear region (see Table |4]l. 
This approach is fully consistent with the first approach explained 
above, which also yields a compact source with = 0.2" since 
FWHM ~ 2Re. 

The difference between the near-IR and optical parameter val- 
ues can be explained mainly by the fact that the extinction is more 
severe in the optical than in the near-IR. Also, in spiral galaxies the 
bulge is mainly composed of late-type stars, hence yielding a more 
prominent bulge in the near-IR than in the optical. Photometry of 
the extracted nuclear source yields Mb ~ —12.24 mag. 



5 DISCUSSION 

5.1 Spectral energy distribution 

NGC 4438 is classified as a LINER 1.9 or a "dwarf" Seyfert 2 
galaxy on the basis of broad Ha emission (Ho et al. I997| |Ken-| 
|ney & Yale|2002| l. We constructed the NGC 4438 SED of the non- 

We did not model the R and / band images because the few pixels cov- 
ering the nucleus were severely saturated. 
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Figure 10. ISAAC Ay grayscale map superposed on TIMM2 NB 1 1.9 /^m 
emission map (contours). Contours are -0.27, 0.24, 0.40, 0.52, 0.61, 0.69, 
0.77, 0.83 and 0.89 times 385 MJy sr"'^. The x- and y-axes are east and 
north offsets in arcsecs from the nucleus. The TIMMI2 astrometry is uncer- 
tain to ~ 5". 




10 15 
log iJ [Hz] 



Figure 11. Spectral energy distribution of distinct nuclear source (AGN) 
in NGC 4438. The solid line is a theoretical SED based on radiative trans- 
fer calculations with a central AGN as the heating source. See text for in- 
formation about data and models. The weak high-energy emission from 
the nucleus of the galaxy suggests that the flows in the accretion disc are 
advection-dominated. 



stellar component using data from radio {u — 4.86 x 10^ Hz) to 
hard X-rays {v = 1.7 x 10^* Hz). 

The near-IR data (log(!//Hz) ~ 14) were obtained from the 
surface brightness decomposition presented in Section]?] The AGN 
component photometry is given in Table |4] and it was computed 
integrating the residual image over a circular aperture of radius 4.5" 
centred on the nucleus. 

The thermal infrared (log(i^/Hz) ~ 13.4) photometry of the 
nuclear source was extracted over a 4.5" circular aperture in the 
TIMMI2 NB 1 1.9 fim image (see Fig.[To|. The result was a mid-IR 
flux density of 0.20 ± 0.04 Jy. Instead of a compact nuclear source. 
Fig. [To] shows that the thermal emission is spatially extended over 
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3.5", which at the Virgo distance corresponds to 270 pc. Thus, the 
nucleus and a circumnuclear dusty arc dominate the mid-IR contin- 
uum emission. 

The soft (log(!//Hz) = 17.6) and hard (log(i//Hz) = 18.2) 
X-ray counterparts were obtained from Chandra data reported in 
Table 4 of |Machacek et aL] ( |2004^ , without any absorption cor- 
rection. The radio continumm (log(j//Hz) = 9.3) data were ob- 
tained from Hummel & Saikia'i fTWT) . This measurement was com- 
puted by integrating the extended emission over a rectangle of sides 
2" X 5". It gives 1.5 ± 0.3 mjy. This value corresponds to an upper 
limit since the nucleus was not clearly detected, and it may contain 
emission from the NW and SE outflows. 

Fig.[TT]shows the SED of NGC 4438. The solid Une is a simple 
theoretical SED taken from ,Siebenmorgen et al.| ( |2004[ l, which is 
based on three parameters: a dusty torus with radius 125 pc, an av- 
erage visual extinction of 4 magnitudes and a compact source with 
bolometric luminosity iboi ~ 5.6 x 10^" L0, which does not rep- 
resent NGC4438's bolometric luminosity because of discrepancy 
between with the model in the X-ray region of the spectrum. The 
AGN is modelled as a power law with monochromatic luminosity 



L{v/Uz 



(i//Hz 



in the wavelength range from 10 A to 



2 fim. The dust is composed of carbon and silicate grains with radii 
between 300 and 2400 A, graphites of radius 10 A and small and 
large PAH components. Fig. [TT] also shows a good match between 
both near-IR and visible data with the theoretical SED. However, 
this model is unable to fit the X-ray data. The sub-Eddington X- 
ray emission Lx ~ 10~® LecM (where LEdd ~ 10^^ erg s^^ is 
the Eddington luminosity for a black hole with mass lO'^ M0) in 
both, the soft and hard bands, suggests that the accretion flows are 
advection-dominated ( [Narayan et al.|1996{|Kenney & Yale|2002| 
[Machacek et al.|2004l >. 

5.2 Outflows and energetics 



From UV observations, |Boselli et al.|p005[ l found that the star for- 
mation in the main body of NGC 4438 is very weak, and is mainly 
composed of old stars, without signs of recent starbursts. Our esti- 
mation of the visual extinction toward the nucleus of Av ~ 4 cor- 
responds to an extinction of 3.3 mag at Ha (Cardel li et al.| 1989) 1, 
implying an extinction-corrected SFR of 0.08 Mq yr~ for the nu- 
cleus. On the other hand, the SFRs inferred from the extinction- 
corrected Pa/3 emission is 1 Mq yr~^ and corresponds to the cir- 
cumnuclear region seen as an arc in Fig. |5ji. These two estimates 
correspond to upper limits since we do not know how much of the 
line emission arises from H II regions near the nucleus. 

The velocity of the northwestern outflow can be estimated 
from the spectrum shown in Fig. |6^. The width of the line con- 
tains contributions from both random motions and bulk velocity. 
The FWHM of the [Fe II] line at 3"- 4" from the nucleus, extracted 
from a 2" aperture is measured to be ~350 km while the mean 
[Fe II] velocity 2" from the nucleus (at the bottom of distinc blob) 
is closer to 500 km (see Fig.|6|. In fact, the mean velocity as a 
function of distance from the nucleus decreases from ~500 km 
at 2" to ~350 km s~^at 2.5". This seems consistent with an ex- 
panding bubble model since if the observed expansion velocity at 
the outer edge of the bubble corresponds to a lower line-of-sight 
velocity than at the middle of the bubble, due to projection effects. 
Hence we adopt a speed of 500 km ^ as a reasonable value for the 
ex pansion velocity, which is a faster outflo w than the one assumed 



(measured from the [Sll] AA6716/6731 doublet ratio). This new 
estimate of the outflow velocity increments the kinetic energy by a 
factor of ~3, yielding an injected kinetic energy into the halo of the 
galaxy of > 10^* erg. This is a lower limit since only ionised gas is 
included in the calculation and we know that there is molecular gas 
also present in the NW bubble (see Fig. [Sj)). The large amount of 
kinetic energy carried by the northwestern outflow along with the 
weak nuclear star formation, suggests that an AGN central engine is 
responsible for much of the optical and near-IR line and continuum 
emission, and not a compact starburst ( [Kenney & Yale|2002[ l. 

5.3 NGC 4438 colours 

The near-infrared colours of the nucleus of NGC 4438 mainly have 
contributions from: a bulge of late-type stars, a non-stellar nuclear 
source, re-radiation from hot dust, and reddening (Kotilainen et alj 
[1992). We measured the colours of each surface brightness com- 
ponent in an aperture of 3", in order to compare with other works 
(there was no substantial variation using an aperture of 8"). Normal 
(inactive) spiral galaxies have colours of 0.6 < J — H < 0.9 and 
H ~ Ks < 0.3 ( |Kotilainen et al.|1992| [Fischer et al.|2006 ^. The 
stellar colours of NGC 4438, J - H = 0.7 and H - = -0.10, 
are clearly in the region occupied by normal spiral nuclei. As ex- 
pected for LINERs and Seyfert 1 galaxies (see two-colour diagrams 
in [Kotilai nen 1993, Forbes et al.||1992} the non-stellar colour of 
NGC 4438 is located in between the regions of AGN and inac- 
tive galaxies, shifting the colours according to a reddened vector 
of Ay ~ 3 — 4. Therefore, the near-IR stellar colours do not show 
evidence for starburst, in agreement with the UV observation of 
[Boselli et al.|P005l l. This seems to be inconsistent with the rela- 
tively high SFR estimated from Pa/3 emission for the circumnuclear 
region (see Section [3.2.3[ ), which would imply that this emission is 
not due to star formation but to the densest gas near the nucleus 
being photoionised by the central AGN. 

5.4 Black hole mass 



[Ferrarese et al.|p006[ > pointed out that a more or less a constant 
fraction of a galaxy bulge mass ends up as a central massive object, 
either a stellar nucleus or a supermassive black hole. The evidence 
discussed in Section [33[ could suggest the presence of a black hole, 
although the radius of its sphere of influence Rbh ~ O.I" is con- 



by [^nney^^kj ((2002| of 300 km s~\ [Kenney & YaIeH2002| 



calculated the mass of ionised gas from the Ha luminosity, which 
corresponds to 3.5 x 10* M© for an electron density of 420 cm"'' 



siderably smaller than our 0.74" resolution. Ferrarese & Merritt 
(2000) found that there is a tight correlation between one of the 
fundamental properties of a galactic bulge, its velocity dispersion 
a, with the mass of its supermassive black holes. Independently, the 
same correlation was found by lGebhardteFallpOOO} , in the same 
year. For NGC 4438 we have calculated the luminosity weighted 
velocity dispersion to be tj = 142 ± 8 km in the central 3 arc- 
sec over the minor-axis, by using the kinematical data (see Fig.[8| 
and by assuming an axial ratio q = 0.57 (see Table[3}. The A^bh-o" 
relation using this value yields a mass for the central black hole of 

log(Mbh/MQ) = 7. 5 ±1.7. 

Recently, Grah am & Driver| ( |2007^ suggested that there is a 
fundamental correlation between the Sersic index n, which is a 
measure of the concentration within the bulge, and the black hole 
mass. The physical interpretation is that steeper and more concen- 
trated bulges (larger n) host more massive black holes. They found 
a Mbh-n log-quadratic relation using updated black hole masses 
and power-law indexes. This relation was computed from a sample 
of 27 galaxies with black hole mass determination from the correla- 
tion between stellar velocity dispersion and virial bulge mass. The 
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Sersic indexes in |Graham & Driver| ( |2007^ were deteimined by per- 
forming a surface brightness decomposition similar to the one used 
in this work. The Mhh-n relation is as tight as the well known cor- 
relation between the stellar velocity dispersion and the black hole 
mass ( [Graham & Driver|2007| l. 

We found a bulge well-represented by a Sersic index n = 1.7. 
This index implies a black hole with mass log(7V4bh/M0) = 
6.70 ± 0.45 for NGC 4438 (for the estimation in each broadband 
see column (16) of Table [3}. The errors in the black hole mass 
were estimated using the intrinsic uncertainties in the [Graham '&\ 
[Driver[ ( |2007^ relation. We did not estimate a black hole mass from 
the bulge parameters inferred from the B-band image, because the 
bulge is highly obscured by dust. Another estimator of Mhi\ is the 
near-IR luminosity ( [Marconi et al.|2003| l, which provides a A^bh- 
ibuigc relation tighter and less sensitive to extinction than those in 
the optical. Column (15) of Table[3]lists the computed black hole 
masses using the near-IR luminosity. It can be seen that these meth- 
ods yield results in reasonable good agreement. 



6 CONCLUSIONS 

In this paper we have presented a study of the nuclear source and 
central environment of the galaxy NGC 4438, based on the results 
of near-infrared ISAAC VLT/ANTU imaging and spectroscopy. 
The main results of this work are listed below. 

(i) We have found extended [Fe II] emission coincident with the 
NW bubble seen in radio continuum. X-rays and optical emission. 
The morphological correlation between [Fe II] and previous radio 
observations, along with the absence of Pa/3 in our spectra of the 
NW bubble, suggest that shocks (perhaps driven by a radio jet) may 
be an important source of excitation of [Fe II] emission. 

(ii) Based on our newly (upwardly) revised estimated of the ex- 
pansion velocity of the bubble, 500 km s~^, and the nuclear star 
formation rate estimated from its emission in Ha, which could be 
up to 0.08 M© yr~^, corrected by the nuclear extinction, we have 
addressed the question whether the outflow has an AGN or a star- 
burst origin. The large kinetic energy associated with the outflow- 
ing gas, > lO^** erg, along with the weak star formation suggest that 
an AGN is more likely to power the outflow. 

(iii) The H2 emission map shows strong molecular hydrogen 
emission around the nucleus which might indicate the presence of 
a molecular torus. 

(iv) Our position-velocity maps showed that the gas emitting 
H2 is kinematically distinct from the gas radiating [Fell] emis- 
sion. Moreover, the nuclear [Fe II] line (FWHM ~ 415 km s~^) is 
broader than the nuclear H2 emission (FWHM ~ 270 kms"^), im- 
plying that the [Fe II] gas is kinematically more disturbed than the 
H2 gas. Furthermore, the molecular gas seems to be more restricted 
to the galactic plane, while [Fe II] extends to higher latitudes from 
the plane of the galaxy. These results suggest that X-ray heating 
from a central AGN may be responsible for the excitation of H2. 

(v) We have applied a two-dimensional surface brightness de- 
composition to the central 30 arcsec of NGC 4438. The best-fitted 
model consists of three components: a Sersic bulge with power-law 
index n ~ 1.7 and ~ 3", an exponential disc with Rh ~ 13", 
and a compact nuclear source, resolved in HST with R^ ~ 0.2". 
The model was evaluated in a generalised elliptical surface with 
discy isophotes and an inclination angle of 69°. 

(vi) We have constructed the spectral energy distribution of the 
nucleus of NGC 4438. It is in agreement with a theoretical SED of 



an AGN obscured by 4 magnitudes, although the X-rays are un- 
derluminous by 6 orders of magnitude (compared to the Eddington 
luminosity for a 10^ M0 black hole), which can be explained on 
the basis of advection-dominated accretion flows. 

(vii) We have derived CO-band stellar line-of-sight kinemat- 
ics along the minor axis of the galaxy. The line-of-sight velocity 
dispersion showed a strong peak on the centre of the galaxy of 
177±8kms'\ 

(viii) A black hole mass of log(A^bh/M0) ~ 7 for the nuclear 
point source could be inferred from the bulge luminosity, the cen- 
tral velocity dispersion and the Sersic index n. 

(ix) We have carried out near-IR photometry of the nucleus, 
finding a point source unresolved by ISAAC at 0.8" resolution with 
Mk = —18.7 and J — K = 0.69. This detection is further evi- 
dence that the central source is associated with an AGN rather than 
a starburst engine. 

(x) We report the brightnesses and colours of the various com- 
ponents in the central region. The dereddened stellar colors (bulge 
and disk) in the nuclear region are typical of an inactive spiral 
galaxy, and show no evidence for a recent starburst. This is fur- 
ther evidence that significant star formation was not triggered in 
the recent collision with M86. 
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APPENDIX A: DATA REDUCTION 

In this Appendix we outline the procedures used in the data re- 
duction process, which included bias subtraction, flat fielding, cor- 
rection of detector defects, cosmic rays removal, and overlapping 
of the dithered frames. The detector has a jump between the two 
halves of the array, caused by imperfect removing of the zero level 
offset due to variations in the bias level ( |Amico et al.||200T| see 
below). 

Flat fielding and bad pixels correction.- The twilight flat field im- 
ages provided by the observatory were sometimes affected by spu- 
rious large scale gradients (specially in the Ks-hand). This could be 
caused by bad ambient conditions, making non-linear effects im- 
portant (L. Tacconi, private communication). We thus divided the 
flat field by a smooth approximation of the large scale variations. 
The smooth flat field was generated by taking the median of the de- 
tector image along the a;-axis, excluding the outer ~ 100 columns, 
which generated a vector column representative of the large scale 
variation. 

Sky and bias residual subtraction.- This is the most important step 
in the reduction. The sky was obtained immediately after the last 
object image, pointing the telescope towards a blank region of the 
sky. The subtraction of the dark current was one of the most im- 
portant problems to solve, because the dark current is known to be 
unstable in ISAAC ( |Amico et al.|200lT l. The bias level of the detec- 
tor is a function of the detector integration time and its illumination. 
The bias level also varies in time, and is more pronounced where 
the readouts start (i.e., rows 1, 2, and rows 513, 514). Due to this 
variation we did not apply the dark subtraction, which yielded im- 
ages containing noticeable bias residuals. We removed these bias 
artifacts by performing a linear extrapolation along the array. The 
procedure was: 

1 . to isolate regions devoid of extended emission, 

2. then, to smooth the regions with a median filter and collapse 
them along the x-axis to obtain single bias residual columns, 

3. finally, perform a linear interpolation between the collapsed 
columns to the entire detector array, generating a smooth bias im- 
age. 

Dithering.- The shifts in the dithered frames were obtained by 
maximizing the covariance between the images as a function of 
astrometric offset (at a sub pixel level). 

Photometric calibration 

We acquired a dedicated set of relatively faint standards from 
the UKIRT telescope system ( Hawarden et al.|2000 ), specifically 
the stars FS6, FS20 (J, H filters) and FS132 (K^ filter). The three 
standard star observations share a similar dithering pattern as the 
object observations, and they were reduced following the procedure 
explained above. The calibration was based on aperture photome- 
try (aperture radius 4.4"). Table [T] summarises the computed zero 
points. We calibrated the narrowband filters using the correspond- 
ing broadband calibration (as recommended in the ISAAC manual, 
[Amico et al.|200l |>. No corrections for atmospheric extinction were 
required because the data were obtained at similar airmass (~ 1.4). 
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Figure Bl. Best-fitted isophotes model to NGC 4438 with constant ellip- 
ticity e = 0.48, angle a = 69° and c = 1.2. This figure also shows the 
convention used for the pixel coordinates X', X, and the angle a. 



al. 



Foreground Galactic extinction obtained from the |Schlegel et| 
1998| maps of dust emission, amounts to E{B -V) = 0.028. 



Table|T[ summarises the magnitudes of extinction for each band- 
pass, computed using Ry = 3.1 and the extinction laws of ]Cardelli| 
[etaLl ( fT989,) . 

Photometric uncertainties were estimated as = (Trms + 
""caiib- The first terms is the rms noise and the second term takes 
into account the calibration accuracy (10 per cent). 

We verified the measured fluxes in our calibrated images using 
2MASS. The result was a 2 — 4 per cent difference between 2MASS 
photometry and our photometry of the standard stars named above, 
and differences of 12, 14 and 10 per cent for JHKs, respectively 
for NGC 4438. This could be caused by an overestimation of the 
sky level in our reduction procedure. 



APPENDIX B: 
MODEL 



SURFACE BRIGHTNESS PROFILE 



This Appendix gives a detailed description of the surface brightness 
model used to fit the NGC 4438 light describution. Our conventions 
and coordinate systems are shown in Fig. |Bl| In order to obtain an 
accurate fit we constructed an elliptical polar-grid, with the coordi- 
nate axis defined by |Athanassoula et aI.H1990t , as follows. 



(1^ 



|2/c 



\y/<i 



2/c^ 



c/2 



(Bl) 



where q is related to the ellipticity by e = 1 — g. c is the pa- 
rameter which allows us to model the generalised elliptical shape: 
c < 1 gives box-like ellipses (boxiness), while for c > 1 we obtain 
shapes approaching to diamonds (diskiness). In general, galaxies 
are inclined at an angle a, so we need to rotate our x- and y-axis 
using a simple transform; X = R{a)'K', where R{a) is the rota- 
tion matrix (see Fig. |Bl| (. The angle a is defined with respect to the 
image pixel coordinate system, increasing counterclockwise. 

The |Sersic| ( |1968^ light profile is very useful for modelling 
elliptical galaxies and steep spiral bulges. It has the form: 



Ib{r) = le exp 



l/n 



(B2) 



where r is the generalised elliptical radius (defined in Equation |Bl[ i, 
le is the intensity at the effective radius R^ and b„ is defined such 
that r(2n) = 2'y{2n,bn) where F and 7 are the complete and 
incomplete gamma functions, respectively. Analytical expressions 
which approximate the value of 6„ give b„ — 2n — 0.33, for 0.5 < 
n < 10. Half of the total luminosity predicted by the profile comes 
from r < Re. The power-law index n describes the shape of the 
light-profile. To calculate the total flux of the Sersic profile, we 
integrated over a generalised elliptical radius r. It gives: 

F^r) = TT he"- Rl ^ 'y{2n,x)I{q,c), (B3) 
where x = b„(r / Re)^^" . The functionT(g, c) is given by 
I{q,c) — / {sxa (f> cos (l)y~ dcj>. (B4) 



This expression accounts for the generalised-elliptical shape. 
When c is greater than 2 ('disky' isophotes) the integral can be 
written as a beta function of the parameter c (see Equation 8 in 
|Peng et al.|2 00?). 

Since NGC 4438 is a spiral galaxy, an exponential law 
is needed to fit the galactic disc. It has the form /d(r) — 
lo exp {—r/Rh), where the parameter Rh is called the disc scale- 
length, lo is the peak of luminosity, and r is the elliptical radius 
(Equation |BTJ. 

A nuclear source can be added at the centre of the galaxy, rep- 
resenting the AGN component, and it may be either a delta function 
(for a point source) or a Sersic component (for resolved sources). 
Finally, the model takes the form: 



/(r) = (/nuc(r) + It{r) + Id{r)) » PSF, 



(B5) 



where J(r) is the surface brightness at a radius r (defined in Equa- 
tion [bTJ, Jnuc represents the nuclear source. It and Id represent 
the bulge and disc components, respectively. The symbol (g) is the 
convolution operator. 

The set of free parameters are: the nuclear source centroid 
(a;ccn, 3/cen), those related to the elliptical geometry: a, q, i and 
c; the Sersic profile power-law index n, the effective radius R^, 
the effective intensity Je; the exponential disc scale-length Rh and 
central luminosity Jo; the nuclear source parameters: n', R'^ and I'^ 
(Sersic case) or the central flux density (delta function case). 

We found the best fitting solution by performing a minimi- 
sation. The main steps of the two-dimensional fitting are roughly 
summarised as follows: 

1 . Select a sub-image, centred on the nucleus. Correct by extinc- 
tion. 

2. Generate a model image based on the initial (or new) condi- 
tions. 

3. Account for the telescope and atmospheric seeing by convolv- 
ing the model with the characteristic PSF. 

4. Evaluate x^. 

5. Iterate from step 2 until convergence is reached. 

6. Generate output images: residual image, original galaxy im- 
age and model. 
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